Preparation and entanglement purification of qubits through Zeno-like measurements 
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A novel method of purification, purification through Zeno-like measurements [H. Nakazato, T. 
Takazawa, and K. Yuasa, Pfiys. Rev. Lett. 90, 060401 (2003)], is discussed extensively and applied 
to a few simple qubit systems. It is explicitly demonstrated how it works and how it is optimized. As 
possible applications, schemes for initialization of multiple qubits and entanglement purification are 
presented, and their efficiency is investigated in detail. Simplicity and flexibility of the idea allow us 
to apply it to various kinds of settings in quantum information and computation, and would provide 
us with useful and practical methods of state preparation. 
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I. INTRODUCTION 

It is usually not seriously discussed in normal text- 
books on quantum mechanics about how to prepare an 
initial state. It is, however, becoming an important sub- 
ject not only from a view point of foundation of quantum 
mechanics, but also from a practical point of view, since 
we are rushing towards experimental realizations of the 
ideas for quantum information and computation 0, Q • 
Without establishing particular initial states assumed in 
several algorithms, we cannot start any processes of the 
attractive ideas. State preparation is one of the key el- 
ements to quantum information processing ^ and 
there are several theoretical proposals |^ H, p and ex- 
perimental attempts H, S IS H ^ • 

In the ideas for quantum information and computa- 
tion, quantum states with high coherence, especially en- 
tangled states, play significant and essential roles. But 
such "clean" states required for quantum information 
technologies are not easily found in nature, since many 
of them are fragile against environmental perturbations 
and suffer from decoherence. Therefore, there would of- 
ten be a demand for preparing a desired pure state out 
of an arbitrary mixed state. Several schemes have been 
proposed for it, which are called "purificatioii," "distilla- 
tion," "concentration," "extraction," etc. (2t [llL[ia |. 

One of the simplest and easiest ways of state prepara- 
tion is to resort to a projective measurement: a quantum 
system shall be in a pure state |</)) after it is measured 
and confirmed to be in the state |<^). Such a strategy is 
not possible, however, in cases where the desired state 
\(f)) cannot be directly measured or where the relevant 
system is not available after the confirmation. This is 
often the case for entangled states, which are the key re- 
sources to quantum information and computation. This 
is why more elaborate purification protocols are required 
and several schemes of entanglement purification/ prepa- 
ration have been proposed ^ 0| . 
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Recently, a novel mechanism to purify quantum states 
has been found and reported: purification through Zeno- 
like measurements |l3j |. A pure state is extracted in 
a quantum system through a series of repeated mea- 
surements (Zeno-like measurements) on another quan- 
tum system in interaction with the former. Since the 
relevant system to be purified is not directly measured 
in this scheme, it would be suitable for such situations 
mentioned above. In this article, we discuss this scheme 
in detail and explore, on a heuristic basis, its potential 
as a useful and effective method of purification of qubits. 
The examples considered here are quite simple but still 
possess potential and practical applicability. 

This article is organized as follows. First, the basic 
framework of the purification is described in a general 
setting, and the conditions for the purification and its 
optimization are summarized in Sec. ^1 where some de- 
tails which are not discussed in the first report are 
included. It is then demonstrated in Sec. lIIII how it works 
and how it can be made optimal in a simplest example, 
i.e., single-quhit purification, and a generalization to a 
multi-qubit case is considered in Sec. II VI which would af- 
ford us a useful method of initialization of multiple qubits. 
One of the interesting applications of the present scheme 
is entanglement purification, which is discussed in Sec.lVI 
and shown to be actually possible. Concluding remarks 
are given in Sec. IVII with some comments on possible ex- 
tensions and future subjects. Appendices A-E are sup- 
plied in order to demonstrate detailed calculations and 
proofs, that are not described in the text. 



II. FRAMEWORK 

Let us recapitulate the framework of the purification 
reported in X^. We consider two quantum systems X 
and A interacting with each other (Fig. The total 
system X+ A is initially in a mixed state gtot , from which 
we try to extract a pure state in A by controlling X. We 
first perform a measurement on X (the zeroth measure- 
ment) to confirm that it is in a state \<p)x. If it is found in 
the state \(j))x, the state of the total system is projected 
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FIG. 1: We repeat measurements on X and purify A. 



by the projection operator 

O = 10): 



1a 



to yield 



^?tot — > Qtot 



Tr(OetotO) 



® QA, 



(1) 



(2) 



where qa = x{<f>\Qtot\(l))x/ Po is the state of A after this 
zeroth confirmation and Pq = Tr{OgtotO) is the proba- 
bihty for this to happen. We then let the total system 
start to evolve under a total Hamiltonian iJtot and repeat 
the same measurement on X at regular time intervals r. 
After N repetitions of successful confirmations, i.e., af- 
ter X is confirmed to be in the state |(^)x successively N 

times, the state of the total system, g[ll{N), is cast into 
the following form: 

= \^)x{^\(^g^l\N), (3a) 
Q^I\n) = {V,{T))^gA{vUT))^/P^^\N), (3b) 



where V0(t), defined by 



(4) 



is a projected time-evolution operator acting on the 
Hilbert space of A, and P^'^\N) is the normalization fac- 
tor. 



e ^^-^)^etot(e^^-^O)^] 



(5) 



Note that we retain only those events where X is found 
in the state \4>)x at every measurement (including the 
zeroth one); other events, resulting in failure to purify A, 
are discarded. The normalization factor P('^)(iV) multi- 
plied by Po, i.e., p(^)(Af) = P('")(A^)Po, is nothing but 
the probability for the successful events and is the prob- 
ability of obtaining the state given in © . 

For definiteness, let us restrict ourselves on finite- 
dimensional systems throughout this article and consider 
the spectral decomposition of the operator V^{t). Since 
the operator V^(t) is not a Hermitian operator, we should 
set up both right and left eigenvalue equations 



V0(r)|M„)A = A„|u„)a, 

A{Vn\V^{T) = X„A{Vn\- 



(6a) 
(6b) 



The eigenvalues A„ are complex in general and bounded 
as 



0< |A„| <1 



(7) 



(see Appendix Here we assume for simplicity that 
the spectrum of the operator V^{t) is not degenerate. In 
such a case, the eigenvectors are orthogonal to each other 
in the sense 



A{Vm\Un)A = S„ 



(8a) 



and form a complete set in the Hilbert space of system 
A, 



X!l^")A(Wn| = 1a, (8b) 



which readily leads to the spectral decomposition of the 
operator V^{t), 



(9) 



(In the following, we also normalize the right eigenvectors 
as a('"„|u„)a = !■) 

Even in a general situation where the spectrum of the 
operator V^{t) is degenerate, the diagonalization Q is 
possible when and only when all the right eigenvectors 
|u„)a are linearly independent of each other and form a 
complete basis 14]. Otherwise, the spectral decomposi- 
tion is not like lO, but in the "Jordan canonical form" 
|T^. The diagonalizability of the operator V^{t) is, how- 
ever, not an essential assumption as clarified in Appendix 

m 

It is now easy to observe the asymptotic behavior of 
the state of A, g^^\N) in (|3 b|) . Since the eigenvalues A„ 
are bounded like (|7I), each term in the expansion 



(10) 



decays out and a single term dominates asymptotically 
as the number of measurements, N, increases. 



(^^(t))^ Xq\uo)a{vo\ as N increases, 
provided 



(11) 



the largest (in magnitude) eigenvalue Xq is 

unique, discrete and nondegenerate. (12) 

[The word "unique" means that there is only one eigen- 
value that has the maximum modulus and "nondegener- 
ate" means that there is only one right eigenvector (and 
a corresponding left eigenvector) belonging to that max- 
imal (in magnitude) eigenvalue.] Thus, the state of A in 
(j3b|) approaches a pure state ]mo)a, 



g^l\N)^\uo)A{uo\ as N- 



oo. 



(13) 
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This is the purification scheme proposed recently : ex- 
traction of a pure state |uo) a through a series of repeated 
measurements on X. Since we repeat measurements (on 
X) as in the case of the quantum Zeno effect [la. we call 
such measurements "Zeno-like measurements" [l^- The 
final pure state |ito)A is the eigenstate of the projected 
time-evolution operator V^{t) belonging to the largest 
(in magnitude) eigenvalue Aq and depends on the pa- 
rameters T, |0)x, and those in the Hamiltonian Htot- It 
is, however, independent of the initial state Qtot- 

The 

pure state |mo)a is extracted from an arbitrary mixed 
state ptot through the Zeno-like measurements. By tun- 
ing such parameters mentioned above, we have a possi- 
bility of extracting a desired pure state |uo)a- 

The above observation shows that the assumption of 
the diagonalizability in ^ is not essential but condition 
(|12|l . i.e., the existence of the unique, discrete and nonde- 
generate largest (in magnitude) eigenvalue Ao, is crucial 
to the purification. For our purification mechanism to 
work, it is crucial that a single state is extracted and 
this is accomplished when these qualifications, i.e., the 
uniqueness of the largest eigenvalue and the nondegen- 
eracy of the eigenvector, are both met. The diagonaliz- 
ability of V0(t) is not relevant to these conditions and is 
not essential to the purification. This point is clarified in 
Appendix IbI 

Furthermore, note the asymptotic behavior of the suc- 
cess probability P^'^^N): it decays asymptotically as 

pM(iV)^ |AoP^PoA(«o|^^A|«o)A 

= |AoP^XA(0i'o|£'tot|'/'Wo)xA as N increases, 

(14) 

where |0uo)xA stands for \(j>)x <Xi |uo)a and xA{(f>U()\ — 
x{4>\ ^ a{uo\- The decay is governed by the eigenvalue 
Ao, and therefore, an efficient purification is possible if 
Ao satisfies the condition 



|Ao|-l, 



(15) 



which suppresses the decay in 1)14(1 to give the final (non- 
vanishing) success probability 



p(^)(iV) 



XA 



{4'Vo\gtot\(t>VQ] 



XA- 



(16) 



It is worth stressing that the condition ((15|l allows us 
to repeat the measurement as many times as we wish 
without running the risk of losing the success probabil- 
ity P'-'^\N). In other words, high fidelity to the target 
state and nonvanishing success probability do not con- 
tradict each other in this scheme, but rather they can be 
achieved simultaneously. At the same time, if the other 
eigenvalues are much smaller than Aq in magnitude. 



purification is achieved quickly. Equations (|15|l and H17|l 
are the conditions for the optimal purification, which we 
try to accomplish by adjusting parameters r, \(f>)x, and 
those in the Hamiltonian i?tot- 

In the following sections, we discuss the above purifi- 
cation scheme in more detail addressing a few specific 
examples, which are so simple but still possess potential 
and practical applications in quantum information and 
computation. 



III. SINGLE-QUBIT PURIFICATION 

Let us first observe how the above mechanism works 
in the simplest example: we consider two qubits (two 
two- level systems) X and A interacting with each other, 
whose total Hamiltonian is given by 



1 



1 



+ gia^a^ + a'^a^), (18) 



where cr^ (z ~ 1,2,3) are the Pauli operators, a± — 
((Ti ±io'2)/2 are the ladder operators, and the frequencies 
rix(A) and the coupling constant g 0) are real param- 
eters. We repeatedly confirm the state of X and purify 
qubit A, i.e., we discuss a purification of a single qubit. 

The four eigenvalues of the total Hamiltonian H^ot in 
((TS|l are given by 



= 0, 



E^^^ ^ (fix + r!A)/2 ± S, 



(19a) 
(19b) 
(19c) 



and the corresponding eigenstates are 

\E^"'>hA = Ui)xA, (20a) 

-ITDxA 



V2\ 



26 



±\/1T^^^^T7^UT)xa), (20b) 



26 

\E^^^)XA = |TT)XA, 

where 

6^^{nx-nA)y4 + g^ 



(20c) 



(21) 



|A„/Ao| < 1 for n^O, 



(17) 



e{g) is the sign function, and |t(i)) is the eigenstate of 
the operator belonging to the eigenvalue -|-1 (—1) with 
the phase convention |t) = cr+|J,). Hence, when the state 
of X, \(j))x, is confirmed repeatedly at time intervals r, 
the relevant operator to be investigated, the projected 
time-evolution operator V0(t), reads 



4 



-iHtatT\ 



V0(t) = x(</'|e' iwx 



sin St sin — 



u>Aai 



cos 6t 



26 



^A . r \ 2'^ 

sm Or cos - 



(22) 



where the state 



is parameterized as 



)x = e-'^^/2 cos -|t)x + e''^/2 sin -|i)x 



(23) 



and the set of angles {6, ip) characterizes the "direction 
of 'spin' X." 

If one of the two eigenvalues of the operator (|22|l is 
larger in magnitude than the other, the condition for pu- 
rification (|12|l is fulfilled, and qubit A is purified into the 
eigenstate |uo)a belonging to the larger (in magnitude) 
eigenvalue Aq. Furthermore, if condition (|15|l . |Ao| = 1, 
is satisfied, we can purify with a nonvanishing success 
probability P^'^^N) xA{4'Vo\gtot\(f>VQ)xA, and another 
condition l|17(l . |Ai/Ao| <C 1, enables us to accomplish 
quick purification. We try to achieve these conditions by 
tuning the parameters. 

The first adjustment for the optimal purification is 



9 = or TT, i.e.. 



|t)x or |i)x (24) 



(see Appendix [n]). Actually, if we choose \<j))x = |T)x, 
the eigenvalues of the projected time-evolution operator 
V^{t) are given by 



Ao 



-i(nx+nA)T 



T — t smOT 

26 



and the eigenvectors belonging to them are 



\Uo)A 

\ui)a 



IT)a, 
IDa, 



a{vo\ = a(TI, 

a(wi| = A(i|- 



(25a) 



(25b) 



It is clear that the magnitude of the eigenvalue Aq is unity 
and that of Ai, 



|Ai| = ^l-(f) sm'6r, 

is less than unity provided 

St ^ nn (n = 1, 2, . . .). 



(26) 



(27) 



Both conditions p2|) and (|15|l are thus satisfied, and ac- 
cording to the theory presented in Sec. ^ we have an 



optimal purification 

'g^l\N)^\]-)A{'i\, 
P(-){N) ^xAiniQiotimxA, 



(28) 



After the repeated confirmations of the state |t)x, qubit 
A is purified into |t)A with a nonvanishing probabil- 
ity XAiTTIftotlTDxA- Similarly, another choice in 
i.e., a series of repeated confirmations of the state 
|J,)x, drives A into |J,)a with a nonvanishing probabil- 
ity XA(ii|£'tot|ii)xA: 



e^A\N)^\l)A{l\, 

P(-\N) xA{ll\gtot\U)xA, 



(TV ^ oo). (29) 



The final success probabihty xA(TT|etot|TT)xA for the for- 
mer choice \(j))x = |T)x or xA(ii|6'tot|ii)xA for the latter 
\(l})x = |i)x means that the target state |TT)xA or |ii)xA 
contained in the initial state gtot is fully extracted. In 
this sense, the purification is optimal. 

The second adjustment is for the fastest purification, 
which is realized by the condition 



(5r (n + l/2)7r (n = 0, 1, . . .), 



(30) 



at which |Ai| in H26|l is the smallest: |Ai| — | fix — f^A 1/2(5. 
We can achieve it by tuning the time interval r, for in- 
stance. 

To be more explicit, let us demonstrate the extraction 
of the pure state |t)A from the initial mixed state 



1 



ftot = |T>x(TI®x(lt>A(TI + U>A(i|). 



(31) 



After X is confirmed to be in the state |t)x successfully 
N times at time intervals r, the state of qubit A and the 
probability for the successful confirmations read 



^sin^^rl^ 



li>A(i| 



J^hM) = lt)A(TI + [l-(g/'5) 

' 1 + [1 ~ {g/ 6)^ Sin' 6t] 

p(^){N)^^{l+[l-{g/6)'sm'6T]''}, 

(32) 

respectively, which clearly confirm the limits H28|) unless 
St = riTT (n = 1,2,...), and the convergences are the 
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fastest when the condition (|30|) is satisfied. (Note that 
XA(TT|f?tot|TT)xA — 1/2 for the initial state considered 
here.) 

In Fig. I^a), the success probabihty P'^'^^N) and the 
so-caUed fidehty to the target state |'Uo)a — |T)a, defined 
by 



F^^HN)^A{u„\g^^\N)\uo)A 
= A(TI^>i^^(iV)IT)A, 



(33) 



are shown as functions of the number of measurements, 
N, for the initial state (|31|l . with the parameters fix = 5, 
^A = 6, g = l,r = 7r/25~ 1.40. Since the condi- 
tion ifT^ . |Ao| = 1, is fulfilled, the decay of the success 
probability P^'^'>{N) is suppressed to yield the finite value 
XA(TT|f?tot|TT)xA — 1/2, and since the time interval r is 
tuned so as to satisfy the condition for the fastest purifi- 



l.Oo 




l.Oo 




N 



FIG. 2: Fidelity F''^\N) and success probability P'-^^N) for 
single-qubit purification. The pure state |t(J,))A is extracted 
from the initial mixed state gtot = |T(J.))x(T(J.)l ^ (IT)a{TI + 
|J.)a(J.|)/2 after repeated confirmations of the state |T(J.))x. 
Parameters are Qx ~ 5, — 6, t — 71/25 ~ 1.40 for (a) 
and fix ~ f^A, T = it/25 ~ 1.57 for (b), in the unit such 
that (7=1. The time interval r is tuned so as to satisfy the 
condition for the fastest purification II3UII in each case. 



cation (|Sn|l (|Ai| ~ 0.45), the pure state |t)A is extracted 
after only iV = 4 or 5 measurements. In an extreme case 
where |Ai| = is possible, the extraction is achieved just 
after one measurement. Such a situation is depicted in 
Fig. IHb) for the same initial state as in Fig. IHa) with 
the parameter set fix = ^Ai 9 — t = — 1.57. 



IV. INITIALIZATION OF MULTIPLE QUBITS 

The single-qubit purification in the previous section is 
too simple but is easily extended for multi-qubit cases. 
In the above example, one may realize that the state 
|TT)xA is an eigenstate of the total Hamiltonian H18|l [see 
(|20|l ] and this is why the optimization condition H15|l . 
|Ao| = 1, is achieved with \(j))x = |T)x and |uo)a = |T)a 
irrespectively of the choice of the time interval r. (The 
same argument apphcs to the case |0)x — \l)x there.) 
In the case of a multi-qubit system X -f A -|- B -|- • • • in 
Fig. 121 with nearest-neighbor interactions. 



1 



n, 



1 + 4 



l + af 



+ 



2 ^2 2 

+ 5xa(^^^^ + 'y^T^) + gAB{<yU^ + a^^r?) 

+ ■■■ (34) 

{gxA,9AB, ... 7^ 0), the state |ttT- • •)xab. . . is an eigen- 
state of this total Hamiltonian i/tot, and it is readily 
expected that the pure state jit. . .)ab. . . is extracted by 
repeated projections onto the state |t)xj with the optimal 
success probability. Similarly, repeated projections onto 
IDx set every qubit into \i) state, i.e., into . .)ab. .. , 
optimally. This would be useful for initialization of mul- 
tiple qubits in a quantum computer. 

In order to make this idea more concrete, let us dis- 
cuss in detail with a three-qubit system X -I- A -f B. The 
important point is whether the condition for the purifi- 
cation (|12|l is achievable, i.e., whether all the eigenvalues 
A„ except for the relevant one Aq, associated with the 
eigenstate |tT)AB (or |ii)AB), can actually be less than 
unity in magnitude. 

For simplicity, we consider the case where ilx = f^A = 
r^B — ^- The eight eigenvalues of the total Hamiltonian 
Htot are given by 





-0, 










p(2) 


= 2fl, 






X 



A 



n ± V2g, 



B 



(35a) 
(35b) 
(35c) 



FIG. 3: A multi-qubit system with nearest-neighbor interac- 
tions. 
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£:(3) = 3n, 

and the corresponding eigenstates are 



XAB 
|4'^>XAB 

|4'^)XAB 



liii) 



XAB, 



cosxUiT>XAB - sinx|Tii)xAB, 

1 



V2 



(sinxUiT)^ 



XAB 



cosxiTii) 



XAB 



± UTi>XAB), 



\E^^^)xAB = cosx|TTi)xAB - sinxUTT) 



XAB, 



(2) 



/XAB 



where 



■^(sinxlTTi)xAB 

± ITiT)xAB), 
ITTT)XAB, 



cosxUTT) 



XAB 



cosx 



V2g 

9XA 



9xA + ffABi 



v5xa+5ab 



smx 



5AB 



V 9xA + 9ab 



(35d) 

(36a) 
(36b) 

i 

(36c) 
(36d) 

i 

(36e) 
(36f) 

(37) 
(38) 



Aiming at initiaUzing qubits A and B into |J,|)ab, we 
repeatedly project X onto the state ||)x at time intervals 
T, and the relevant operator to be investigated reads 

F4(r)^x(i|e-^^-ni>x 

= Ui)AB(ii| 

+ |Ti>AB(n|e-''"" cosv^gr 

+ UT)AB(iT|e-^''"(cos2x + sin^xcos V2gT) 

^ t\n)AB{.me-'"^ smxsmV2gT 

-«UT)AB(Ti|e"*''"sinxsin\/2.gT 

+ |TT)AB(TT|e-'''""(sin2x + cos^xcos V2gT). 

(39) 

The target state HDab is an eigenstate of this operator 
belonging to the eigenvalue = 1, which satisfies the 
optimization condition (|f 5(1 . and the other three eigen- 
values are give by 



Ah 



e-f^-fcos^^-sin^xsin^ 



V2 



V2 



EL 

V2 



■ a • / A ■ ? gT , ■ 2 n 9T 

=F sm \ I cos^x sm — 4 sm x cos^ 



An=e-2'^^Ml-2cos^Xsm2^ 



V2j' 
(40a) 

(40b) 



If these three eigenvalues are all less than unity in mag- 
nitude, the condition for the purification l|12|) is satisfied, 
and the initialized state HDab is extracted from an arbi- 
trary mixed state gtot , with a nonvanishing success prob- 
ability p(")(iV) ^ xAB(iii|^'tot|iii)xAB. (Note that the 



-H 




FIG. 4: Magnitudes of the eigenvalues A+ (solid line), A_ 
(dashed line), and A|| (dotted line) in 14011 . as functions of 
V^gr. In this figure, we set grxA = Qab- Note that |A+| = |A_ 
within each range 2nn - C < V2gT < 2n-K + ( (n = 0, 1, . . .), 
where ^ (0 < C < ■"") is defined by tan((^/2) = 2| sinxl/ cos^x, 
and ( ~ 0.787r when gxA ~ 3ab. 
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FIG. 5: Fidelity F^^\N) and success probability P'^'(iV) 
for two-qubit initialization. Through the repeated confirma- 
tions of the state |i)x, qubits A and B are initialized into 
i.i)AB from the thermal equilibrium state of the total system 
at temperature T, i.e., ^itot oc e"^^'"* with /3 = (ksT)'^ . 
Parameters are gxA ~ <?ab = f, f2 = 2, r = ~ 1.73, 

ksT = f3~^ = I. The time interval r is tuned so as to make 
max(|A+|, |A-|, |A|||) the smallest, which is for the fastest ini- 
tialization (see Fig.^J- 



left eigenvector belonging to the eigenvalue Xn is AB(iil-) 
Such a situation is realized provided 



V2gT^mT (n = l,2,...), 



(41) 



which is clearly seen from Fig. ^ and a proof in Ap- 
pendix ^ The final success probability P<^''\N) 
XABdiilftotliiDxAB is again optimal, in the sense that 
the target state |iii)xAB contained in the initial state 
gtot is fully extracted. 



7 




FIG. 6: We repeat measurements on qubit X and extract one 
of the Bell states, |*~)ab = (ITDab - |iT)AB)/V2, in qubits 
A and B. 



The above argument reveals the possibihty of initial- 
ization at least for two qubits. Initialization of two qubits 
into |||)ab from the thermal equilibrium state of the to- 
tal system at temperature T, i.e., gtot oc e~''^'°' with 
(3 = (fcer)"^, is demonstrated in Fig. O Note that it is 
effective when > V^g, since in such a case, |iii)xAB 
is the ground state of the total system. The analytic 
formula for the final success probability is P^'^^(oo) — 
[1 + (e-^^ + e-2/3^2)(i + 2cosh 72/35) + e-^'^^]~\ 

It is natural to expect that the same mechanism also 
works for systems with more qubits as in Fig. |31 It is 
hard to imagine that the magnitudes of eigenvalues of 
Vi{t) other than the relevant one Xn... (whose magni- 
tude is unity) is also unity irrespectively of the values of 
parameters. Further detailed investigations on its effi- 
ciency, robustness, and so on, will certainly clarify the 
possibility of a new useful procedure for initializing mul- 
tiple qubits. 



ENTANGLEMENT PURIFICATION 



tanglement, and therefore, it is interesting and important 
to examine whether the present scheme can realize entan- 
glement purification/preparation. We show, in this sec- 
tion, that it is actually possible. In order to demonstrate 
it explicitly, let us discuss a simple Hamiltonian 



-fftnt — fi- 



n- 



,X A 



n- 



(42) 



The control qubit X is coupled to qubits A and B as in 
Fig. El We confirm X to be in the state \4>)yi repeatedly 
at time intervals t and end up with an extraction of an 
entanglement between A and B, which are initially in a 
mixed state ptot- 

The spectrum of the total Hamiltonian i?tot is already 
given in H35|l with g replaced by and the eigenstates 
are 



\E, 



.(1) 



XAB — 



'0 /XAB 



\E. 



± /XAB — 



1 



i)xAB, 
XAB, 



-^[Uvl/+)xAB±e(5)|ni)xAB] 



|4'^>XAB = 



i(2)\ 



= IT*" 



/XAB, 



f ';XAB = -^[|T«'+)XAB ± e(5)UTT)xAB] 
£;^'^)XAB - |TTT>XAB. 



(43a) 
(43b) 

(43c) 

(43d) 
(43e) 

(43f) 



One of the most significant issues in the field of quan- 
tum information and computation is how to prepare en- 



The relevant projected time-evolution operator is given, 
in this case, by 



y4r)EEx(0|e-*^-10)x 

= |vI/-)AB(*-|e-"- (^sin^^ +e-"^cos20 

+ Ui>AB(ii| (^sin2^+e-*f^-cosV2grcos2 

+ |*+)AB(vI'+|e-*"^cosV2gT (^sin^^ + g-^^^ cos^^ 

+ |TT)AB(TT|e-''''" (^cos\/25rsin2^ + e-'"^cos2 

- i (Ui)AB(*+|e"^ + |*+)AB(ii|e-'^) e-'"^ sin V2.gr sin ^ cos ^ 

- I (|tT)AB(*+|e-*'^ + |*+)AB(TT|e'^) e-2^f^^ sin V^gr sin ^ cos ^, (44) 
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where |^'^)ab are the two of the four Bell states 
|*±)ab = (ITDab ± |iT)AB)/V2, |$±)ab = (ITT)ab ± 
|ii)AB)/V2, and |(?!))x is parameterized as in (j2Sll- Since 
the Hamiltonian (|42|l is symmetric under the exchange 
between A and B, V^(r) splits into two sectors: the sin- 
glet sector and the triplet one. The singlet state |^'~)ab 
is apparently one of the four eigenstates of (r) belong- 
ing to the eigenvalue 



A^- = e 



sm - 
2 



cos^- 
2 



(45) 



and hence, we can extract an entangled state, i.e., the 
Bell state |\1/~)ab, after a number of measurements on 
X, provided (i) the eigenvalue is larger in magni- 
tude than any other eigenvalues. Furthermore, if (ii) 
condition (fT^ . i.e., |A^-| = 1, is achieved, |^'^)ab 
is extracted with an optimal probability P'-'^\N) 
XAB('^^~|f?tot|'^^~)xAB, which is again optimal in the 
same sense as in the preceding examples, i.e., the target 
entangled state |4'~)ab contained in the initial state ptot 
has been fully extracted. 

Requirement (ii) is fulfilled by the choice of the pa- 
rameters as |r2|r — 2mT (n = 0, 1, . . .) or sm9 = 0, but 
the latter choice violates requirement (i). It is, therefore, 
necessary that 

|r!|T = 2n^ (n-0,1,...) and |0)x |T)x, |i)x. 

(46a) 

(Note that the first condition is automatically satisfied 
without tuning the time interval r, when Q ~ 0.) Fur- 
thermore, one can prove as in Appendix IeI that require- 
ment (i) is met, under the condition 14()all . provided 



(46b) 



\g\T/V2^m7:/2 (m=l,2,...). 



The existence of such a parameter set satisfying (|46ll ex- 
plicitly discloses the possibility of extracting entangle- 
ment through Zeno-like measurements. 
In the case of the choice 

l'/'>x = H)x^^(|T)x + U)x), (47) 

for example, the four eigenvalues are given by 
9T 



A 



1, 



A<£- 



V2' 



(48a) 



A+ = 1 — i sin (3 sin ± e(g), / 1 
2 ^2 V V2 ^' 



9cos2^ ) , 
(48b) 



whose magnitudes behave as in Fig. ^ but with A|| re- 
placed by A$- . 

The extraction of the entangled state |\E'~)ab is 
demonstrated in Fig. [TJ from a product state gtot — 
|— >)x(^|'8) |T)a(TI® li)B(i| and from the thermal equilib- 
rium state ptot oc e~^'^*-°*- at temperature T = {k^(3)~^ . 





FIG. 7: Fidelity F^^\N) and success probability P'-^^N) for 
entanglement purification. The entangled state |^~)ab is ex- 
tracted from (a) a product state ^itot = |— '•)x(— »| <8) |T)a(TI ® 
|J.)B(i.i and (b) the thermal state Qtot oc e~^^''°* at temper- 
ature T = (fcB/3)~^, through repeated confirmations of the 
state l^)x. Parameters are — 0, t — O.Svr ~ 1.57 for (a), 
and n = 0, r = C/\/2 ~ 1.73, fceT = /^'^ = 00 for (b), in 
the unit such that g — 1, where ^ is defined in the caption 
of Fig. |1| For the initial thermal state in (b) with f2 = 0, 
the success probability for the zeroth confirmation is given by 
P'^'(O) — 1/2 for any set of parameters {9, ip, g,T), and the 
final value P''^'(oo) = [8 cosh^ {13 g/\^)]~^ becomes largest at 
kBT/\g\^{P\g\)-' ^<x>. 



VI. CONCLUDING REMARKS 

The examples presented in this article demonstrate 
how the present purification scheme works, and suggest a 
few potential applications, even though the analyses are 
heuristically based and no general "optimization" theory 
or strategy has been given. Remarkable features of the 
scheme are summarized as follows, (i) The first point 
is the simplicity. Many of the other proposed proce- 
dures are composed of several steps with different op- 
erations, such as rotation, CNOT operation, and mea- 
surement 0, In the present scheme, on the other 
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hand, one has only to repeat one and the same mea- 
surement, (ii) Furthermore, the "optimal" success prob- 
ability is possible in the sense that the target state con- 
tained in the initial state is fully extracted. In several 
other methods [2, 0| , on the contrary, it decays to zero 
as the fidelity approaches unity |18| |. (iii) The number 
of measurements required for purification is considerably 
reduced by appropriate choices of parameters, and pu- 
rification is attainable after only a few steps. 

Another point to be stressed is the flexibility. While 
many of the other schemes 2, 11] arc designed for specific 
systems, the framework is presented in Sec. ^] on a gen- 
eral setting, and there are diverse systems and purposes 
which fit the present scheme. We have already observed, 
in this article, two different applications on the same idea: 
initialization and entanglement purification. Additional 
ideas or slight modifications to the basic scheme would 
provide us with various methods of state preparation. 
An interesting extension of the present scheme is for ex- 
traction of entanglement between two spatially- separated 
qubits 0,ll[illi2, which is often necessary for quantum 
communication, quantum teleportation, and so on. (The 
original protocols for entanglement purification 0, ^ 
are aimed at this purpose.) It is actually possible and 
will be reported elsewhere [T3 . One of the other possible 
extensions is to go beyond a method of extracting quan- 
tum state. It would be interesting, for example, if we 
could find a novel method of transferring quantum state 
[2(tI I21II rather than extracting it. 

In this article, only qubit systems, i.e., finite-dimen- 
sional systems, have been discussed. One has to keep 
in mind that the condition H12|l plays a crucial role in 
the present purification scheme. If this condition is met, 
however, it works for infinite-dimensional ones as well. 
In fact, a harmonic oscillator, which has an infinite num- 
ber of energy levels, can be purified through the present 
method, which is explicitly demonstrated in This 
also shows the broad range of applicability of the scheme. 
It is not obvious, however, whether one can purify sys- 
tems with continuous spectra, since they seem, at first 
sight, unlikely to satisfy the condition for purification 
(|12|l . especially the discreteness of the eigenvalues. This 
point is one of the interesting future subjects, since it 
would be required in some cases to purify quantum states 
in the presence of environmental systems, namely, under 
dissipation and/or dephasing. 

The simplicity and the efficiency mentioned above 
would facilitate practical experimental applications of the 
present scheme. The flexibility allows one to apply it to 
various kinds of systems intended for quantum informa- 
tion and computation, such as optical setups |l2j |. ion- 
trap systems 0, |^, solid-state quantum computers l^, 
and so on. In practice, one should face many unwanted 
factors, and robustness of the method against them is 
crucial. In the present scheme, it is often required to 
tune certain parameters in order to extract a desired pure 
state, and it is an important subject to clarify how pre- 
cise the tuning should be and how much error the method 



suffers from when the parameters are mistuned. It is also 
a remained issue to explore how ideal projective measure- 
ments are realized in actual experiments. Investigations 
on these points are now in progress. 
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APPENDIX A: BOUND ON THE EIGENVALUES 

OF V4,{r) 

Let us prove that the eigenvalues X„ of the projected 
time-evolution operator V^{t) are bounded as in l|7|l. 
For an arbitrary state of A, say |V')a, 

o<\mTmAf 

< ||e-'^--|0)x|^>A|P 

= 1- (Al) 

Hence, by setting |^)a — \u„)a [a right eigenvector of 
the operator V^{t)] and noting || V0(r)|u„)A|P = |A„p, 
we obtain the inequality Q. As is clear from this proof, 
the bound Q reflects unitarity of the time-evolution op- 
erator e"*^'"'"^. 



APPENDIX B: A NONDIAGONALIZABLE V^(r) 
CASE 

It is assumed in Sec. ^ that the projected time- 
evolution operator V^{t) is diagonalized like but it 
is not the case if some of its eigenvalues are degenerated. 
Here we show, however, that the assumption of the diag- 
onalizability is not essential to the purification. 

When an eigenvalue A„ of the (finite-dimensional) op- 
erator V(I,{t) is M„-fold degenerate, there do not always 
exist Mn linearly independent eigenvectors. This fact 
spoils the diagonalizability of the operator V^{t). There 
exist dn (< Mn) linearly independent right eigenvectors 

\uii'')a {k = l,...,d„) belonging to the eigenvalue A„ 
{dn is called "dimension of the eigenspace" ) , and one can 
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find Mn — dn linearly independent "generalized eigenvec- 
tors" \uI!'^)a (fc = (i„ + 1, . . . , Mn) which are subjected 
to the conditions 

y4r)|«W)A-A„|«W)A + k'^-i))A 

(fc = d„ + l,...,M„) (Bl) 

and linearly independent of the eigenvectors |m1'^'')a {k = 
1, . . . , d„) [l|. The right vectors (fc = 1, . . . , M„) 

then form a complete set within the subspace associated 
with the eigenvalue A„, and there exist corresponding 
left vectors A(wi'^''| {k ~ 1,...,M„), which satisfy the 
orthonormality 



where 



and completeness conditions 



fc=l 



The operator V0(t) is now expanded as 



with 



fc=d„ + l 



(B2a) 



(B2b) 



(B3a) 



(B3b) 



which is the most general form of spectral decomposition 
and is called "Jordan canonical form" Note the 

relations 



and 



VrnVn^O {m ^ n) , 



(B4a) 
(B4b) 
(B4c) 



(B4d) 



From the spectral decomposition (jB3|) . it is easily de- 
duced that 

, min(N,M„-dn) , 



r=l 



M„ 

E I 

k—dji-\-r 



(k-r) 



)a{v.. 



(fc)| 



(B5b) 



Therefore, if the largest (in magnitude) eigenvalue is 
unique, which is denoted by Aq, and nondegenerate (i.e., 
Mq = 1, do = 1, Dq — 0), the single term in the ex- 
pansion ljB5a|l again dominates asymptotically like Hll() 
(note that nCt ^ N'^ /r\ for large N), which leads to 
the same conclusion as (I13|l . The purification does not 
suffer from degeneracy in the other eigenvalues than the 
largest (in magnitude) one Aq. The crucial condition to 
the purification is 

APPENDIX C: OPTIMIZATION OF THE 
SINGLE-QUBIT PURIFICATION 



We show here that the condition (|24|l together with 
()27|l is the necessary and sufficient condition for the op- 
timal purification with both (|12|) and (|15|l for model l|18|) . 

First, we try to achieve the upper bound in the in- 
equality lPT|l . i.e., |1V0(t)|?/')a|| = 1, in model lO. If 
such a state |V')a is found and is an eigenstate of the op- 
erator V0(r), say |m„)a, we have |A„| — 1. As is easily 
seen from (|A1|) . the equality holds only when 



-iHtotT\ 



:|V')a = 



(CI) 



is satisfied, where |(/>j.)x is a vector perpendicular to |(/))x 
in 1231), i.e.. 



|</>±)x = e-^^/^ sin ^|T)x " cos -U)x. (C2) 



Equation l|Cl|l means that the operator V^{t) = x{4>± \ x 
e~*^'°''^|0)x should have a zero eigenvalue, and hence 



dety/(T) =0. 



(C3) 



(B5a) For model (fH^ . the operator Vi'-(r) reads 



_ i(Jix+JU)r/2 



cos St + i 



. f^X — f^A 



26 



sin Si 



sm — cos — 
2 2 
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-u)Aai 



26 



sm - cos - 
2 2 



2 2 



and 



detvHT) = -]e 



26 



sin (5t 



+ (|) sin^^r 



(C5) 



r 



Condition HC1|I . namely (jC3p . is hence reduced to 

sin 61 = (C6a) 



or 



cos (5t = ±1 and 



„i(nx+nA)r/2 



±1. 



(C6b) 



In the first case (IC6a|l . both conditions (|12|l and p5|) 
are satisfied unless 5t — nir [n = 1,2,.. .) as is ex- 
plained around l|24|) - (|27|l . In the second case (|C6bp . on 
the other hand, the projected time-evolution operator 
reads V^{t) — 1a and the eigenvalue Aq = 1 is degener- 
ated, i.e., condition \i2\ is not fulfilled. Therefore, the 
necessary and sufiicient condition for the optimal purifi- 
cation in model (|18|l is given by the first choice ljC6a|) 
[i.e., (EH)] with ||2Zl. 



APPENDIX D: CONDITION FOR THE 
TWO-QUBIT INITIALIZATION 

We here outline the proof of the necessary and suffi- 
cient condition for the optimal two-qubit initialization, 
Eq. H41|l . in Sec. IIVI What we have to show is how to 
make the eigenvalues \± and A|| in H4U|) all less than 
unity in magnitude. 

The eigenvalues A± are the solutions to an eigenvalue 
equation 



[Xe^^-f - 2 ( cos2^ - sin^xsin"^ ) (Ae'^'^) 



2 • 2 5''" \ / \ iQ.T\ 
V2) 



1 - 2cos^Ysin2-^ = 0. 

V2 



(Dl) 

We clarify when this equation has a solution whose 
magnitude is unity. Seeking such a solution, we insert 
A = e^'^'^e'® into HDlfl to obtain the conditions 



sin e I cos e - cos^ ^ + sin^x sin^ ^1=0, (D2a) 



2 



V2 



V2 



cos 9 cos 6 - cos ^ + sin'^x sin ^ 



cos^Xsin^^ = 0, (D2b) 
v2 



which are reduced to 

. gT 

sm —pz — 

or 



and cos = 1 



cos = and cos 6 
V2 



-1. 



(D3a) 



(D3b) 



(Note that cosx 7^ and sinx 7^ 0, since it is assumed 
that 5xA,ffAB ^0.) It is easy to see from H40b|) that 
we have |A||| = 1 when sin(gr/V2) = 0, and in sum- 
mary, the magnitude of one of the eigenvalues A± and 
A II becomes unity only when 



^/2gT — rnr (n — 1,2, . . .). 



(D4) 



The condition for the initialization in Sec. IIVI is thus 
proved to be (jUl). See also Fig. El 



APPENDIX E: CONDITION FOR THE 
ENTANGLEMENT PURIFICATION 

The necessary and sufficient condition H46|l for the en- 
tanglement purification in Sec. is proved in a similar 
manner to that in Appendix IdI 

The eigenvalues A± and A$- under the condition H4t)a|) 
are the solutions to an eigenvalue equation 



K — sm 



9T 

V2 



2" • 9''' 
K — COS — sm 
2 



V2 

El 

V2 



. 2^ ■ 9T 

K — sm — sm 
2 



V2 



2 sin^ ^ cos^ ^ cos^ ^ in - ^sin^ 1 =0 



V2) 



(El) 



with K = (1 — A)/[2sin(gT/\/2)]. Seeking a solution A 
with unit magnitude, we insert A = e*® into this equation 
to obtain 



sin Q 



2cos^Q- ( 3-4sin2^ 1 cosO + l 



- ( 2 - - sin^e sin^ 



V2 

9T -2 9T 

^ 2 - sm —= 
V2\ V2 



= 0, 
(E2a) 



cos 9 



2cos^e - ( 3- 4sin2-^ ) cosG 
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— 2 — - sm 6^ sm —= 2 — sm —= 

V 2 ; V2V 72 

_l_2sii/^-isin20sin2:^ f2-3sin2^ 
V2 2 V2 V \/2 



0, 



(E2b) 



or 



or 

cos — = = and cos O 
V2 



(E3b) 



which are reduced to 

9T 



sin -2— = and cos 6 = 1 (E3a) 
v2 



Extraction of entanglement is not possible when (|E3a|) 
or ljE3b|l is satisfied, and therefore, the condition for the 
entanglement pmification in Sec. is given by (j^ . 
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